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The dynamics of an F-center created by an oxygen vacancy on the TiO2(110) rutile surface has 
been investigated using ab initio molecular dynamics. These simulations uncover a truly complex, 
time-dependent behavior of fluctuating electron localization topologies in the vicinity of the oxygen 
vacancy. Although the two excess electrons are found to populate preferentially the second subsur- 
face layer, they occasionally visit surface sites and also the third subsurface layer. This dynamical 
behavior of the excess charge explains hitherto conflicting interpretations of both theoretical findings 
and experimental data. 
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Titanium dioxide (Ti02) is one of the most thoroughly 
investigated metal oxides, due to its broad range of uses 
in several key technologies including heterogeneous catal- 
ysis, pigment materials, photocatalysis, and energy pro- 
duction, to name but a few [J-Q. It is well known that 
bulk and surface defects govern the properties of tita- 
nia, and are thus of fundamental importance in virtu- 
ally all its applications The most common point 
defects on the TiO2(110) rutile surface are oxygen va- 
cancies (O v ) in the two-fold coordinated O rows and Ti 
inter st it ials @, Ej. In particular, removal of an O atom 
gives rise to two excess electrons and the appearance of 
new electronic states in the band gap at about 0.7-0.9 eV 
below the conduction band edge creating an F-center 

Although the two excess electrons can in principle 
be localized on any Ti atom, they are believed to pref- 
erentially occupy specific Ti-3d orbitals, thus formally 
creating Ti 3+ sites [Io|, [l2[. In stark contrast, recent 
experiments [13| suggest a qualitatively different view- 
point: charge localization is found to be more disperse, 
with the excess electrons being shared by several surface 
and subsurface Ti ions. Furthermore, STM/STS experi- 
ments have revealed charge derealization involving more 
than ten Ti sites [14j . 

Unfortunately, different computational methods yield 
conflicting results [ll[. Local/semilocal density function- 
al (LDA/GGA) predict a rather delocalized defect level 
for O vacancies on TiO2(110) with an energy right at 
the bottom of the conduction band [ll[. However, it is 
well known that such functionals bias against localization 
on strongly correlated d-states, and hence alternative 
methodologies are welcome. Recent studies of defective 
Ti02 surfaces [151421] have focused on "pragmatic and 
practical" correction schemes using hybrid functionals or 
a Hubbard correction. Although both schemes yield the 
expected gap states, they each predict vastly different 
localization topologies of the excess charge. 

Using B3LYP on a c(4x2) slab with an O vacancy, the 
defect charge is found to be localized on d-orbitals of two 
surface Ti atoms [15j . In particular, one unpaired elec- 



tron is found on the under-coordinated Ti(ll) site, while 
the other moves to an adjacent five-fold coordinated 
Ti5 C atom, such as Ti(7); see Fig. [T] for our site label- 
ing scheme. By contrast, LDA/GGA + U studies 16l-[21| 
on the reduced TiO2(110) surface have reported charge 
localization on different surface and/or subsurface sites. 
For instance, a combination of surface and subsurface lo- 
calization immediately beneath the defect on Ti(ll) and 
Ti(27) (see Fig. [1]) has been found [16[. However, the 
results are reported to be strongly dependent on the su- 
percell size. Using a (2x1) cell, the electrons are found 
at Ti(7) and Ti(23), while using a (4x1) cell complete 
subsurface localization is observed [2l[ at Ti(23) and 
Ti(39). On the other hand, a (4x2) cell yields local- 
ization of the electrons on Ti(ll) and Ti(12) when using 
U > 4.2 eV, whereas smaller values of U lead instead to 
derealization [19] . Interestingly, some recent GGA + U 
studies 22|, |23[ which focussed on the intrinsic electron 
transport in titania provide useful hints: Electron hop- 
ping in defect-free TiC>2 bulk [22] is described by a po- 
laron localized at a Ti 3+ site, which can hop to an ad- 
jacent Ti 4+ with activation energies of only « 0.1 eV. 
Furthermore, static calculations of ideal TiO2(110) sur- 
faces with a single excess electron suggest that the most 
favorable trapping sites are subsurface ions that lie below 
rows of Ti 5c atoms on the surface . 

In summary, the existing results of electronic struc- 
ture calculations which suggest highly localized excess 
electrons at very specific Ti sites are at odds with the 
current experimental picture of a more "fuzzy" scenario 
involving many sites and delocalized excess charge. How- 
ever, low activation energies (of the order of 0.1 eV) and 
the multitude of topologies for the localization of excess 
charge that have been previously revealed suggest that 
the energy landscape might be relatively flat and thus 
prone to thermal fluctuations which could affect localiza- 
tion scenarios in a dynamical sense. 

In this Letter, we describe extensive GGA + U ab ini- 
tio molecular dynamics (AIMD) simulations [24| of the 
excess charge (de-) localization dynamics as induced by 
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FIG. 1: (a) Ball and stick model of the defective TiO 2 (H0) 
surface. Red and blue spheres are O and Ti atoms, respec- 
tively. Panels (b), (c), and (d) depict the spin density (at 
0.005 e/A 3 ) of three configurations from Table [I] with differ- 
ent charge localization topologies. 



temperature. We find that the defect charge created by 
an O vacancy on TiO2(110) is dynamically sharedby sev- 
eral subsurface and surface Ti sites, with a dominant 
contribution from particular second layer subsurface sites 
which do not belong to Ob rows. Thermal fluctuations 
of the titania lattice allow the excess charge to probe - 
and thus to populate - various local (electronic structure) 
minima of similar energy but vastly different localization 
topologies. Our findings demonstrate the need to go be- 
yond static optimization in order to uncover the dynam- 
ical nature of the phenomenon. These results strongly 
support the proposal of a delocalized nature of the ex- 
cess charge for such defects, and in addition, they also 
reconcile the hitherto contradictory viewpoints from dif- 
ferent static electronic structure calculations. 

The reduced TiO2(110) surfaces have been modeled 
by four 0-Ti202~0 tri-layer (4x2) supercell slabs, sep- 
arated by more than 10 A (see FigH]). The bottom of 
the slabs was passivated with pseudohydrogen atoms of 
nuclear charge +4/3 and +2/3 in order to achieve well 
converged results, which is our "standard setup" [llj. All 
calculations have been performed using spin-polarized 
GGA + U, PBE 0, and ultrasoft pseudopotentials [27j | . 
Extending our previous work [25|, we have implemented 
GGA + U into CPMDS using the self-consistent linear 
response approach [29|, |30[ to compute the Hubbard pa- 
rameter, yielding U = 4.2 eV for our setup; as usual the 
localization details will depend on the particular value 
of U ', but the qualitative features reported in this paper 



TABLE I: Relative energies (in eV) and Ti(ll)-0 bond 
lengths (in A) of reduced TiO2(110) with different localiza- 
tion of the O v -induced excess charge. The last three columns 
show on which layer (s) the two excess electrons are localized. 
Bold labels refer to the configurations in Fig. [TJ 



Ti 



3+ 



Ti 3+ AE Ti(ll)-0 l st -layer 2 nd -layer 3 rd -layer 



Ti(6) Ti(32) 0.23 1.87 1 e 1 e e 

Ti(15) Ti(21) 0.28 1.87 1 e 1 e e 

Ti(ll) Ti(29) 0.77 1.92 1 e 1 e e 

Ti(24) Ti(31) 0.00 1.87 e 2 e e 

Ti(22) Ti(29) 0.04 1.87 e 2 e e 

Ti(24) Ti(30) 0.16 1.86 e 2 e e 

Ti(27) Ti(29) 0.22 1.86 e 2 e e 

Ti(30) Ti(40) 0.36 1.87 e 1 e 1 e 



were checked to be stable upon reasonable variation. The 
occupations of the d-orbitals are calculated using atomic- 
like wavefunction projectors. The lowest tri-layer atoms 
are constrained to their equilibrium positions, while all 
other atoms are free to move. The AIMD simulations 
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were on the order of 10 ps in length, and used the Car- 
Parrinello propagation technique [3l| with a fictitious 
electron mass of 700 a.u. and a time step of 0.145 fs. 

In order to reveal the nature and distribution of the 
defect charge, we performed the PBE + U simulations at 
various temperatures, ranging from 700 to 1000 K. This 
is well below the melting point of almost 2000 K, but 
sufficiently far above the ambient temperature that the 
phonon dynamics is accelerated and thus the sampling 
of the potential energy surface (PES) is enhanced on the 
picosecond AIMD time scale. Two complementary sce- 
narios, characterized by different defect charge localiza- 
tion topologies, were employed in order to provide two 
distinct sets of initial conditions. One scenario repre- 
sents a case where the two excess electrons are trapped 
at two second layer sites, Ti(24) and Ti(30), under Tis c 
rows. The second scenario corresponds to localization on 
an under-coordinated site in the first layer, Ti(ll), and 
on Ti(29) in the second layer under a Ti 5c row. In both 
cases the paramagnetic triplet state is preferred and the 
structure with defect charge localized on the second layer 
is roughly 0.6 eV more stable (see Table H]). 

After equilibrating both setups at 300 K for about 2 ps 
each, the systems were heated to their respective simu- 
lation temperatures. The charge localization dynamics 
was monitored by computing, as a function of time, the 
occupation matrix of each Ti d-a and d-f3 spin orbital 
along the trajectories. We will first consider the effect 
induced by thermal fluctuations on the setup with all 
excess charge initially localized in the second layer. As 
demonstrated by Fig. Efa), the defect charge transfers 
from sites Ti(24) and Ti(31) to Ti(21) and Ti(32) in less 
than 0.5 ps; we note in passing that a qualitatively similar 
scenario is observed using a smaller (2x1) cell. Approx- 
imately 30 such charge migration events were observed 
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FIG. 2: (a) Dynamics of the fractional occupation of partic- 
ular Ti d-orbitals during a time fragment of the PBE + U 
simulation carried out at 1000 K; at t = the charge is lo- 
calized in the second subsurface layer at Ti(24) and Ti(31); 
populations of about 1 and correspond to Ti 3+ and Ti 4+ 
charge states, respectively, (b) Same as (a) but the charge is 
initially localized in the first layer at Ti(ll) and subsurface at 
Ti(29) at 700 K. (c) Distribution function of the average pop- 
ulation of all available Ti-sites by the two excess electrons 
obtained from the full simulation underlying (a); inset pro- 
vides layer-averaged populations, (d) Electronic DOS (a/ (3 
spins: red/black lines) of the three structures from Fig. [T] 
peaks in the band gap correspond to reduced Ti-sites; the 
top of the valence band is set to eV, and the red area filling 
the gap states (scaled by 10 to enhance visibility) stems from 
Ti 3+ sites as labeled. 



during the 6 ps simulation, although only a few of them 
involve charge leaving the second subsurface layer and 
moving to the first or third layer. This implies that the 
defect charge spends most of the time at second-layer Ti 
ions, but translocates very rapidly from a given Ti site to 
an adjacent one within the same subsurface row. We ob- 
served no charge transfer to Ti sites belonging to adjacent 
rows. Other interesting electronic topologies explored by 
the dynamics are cases in which the two excess electrons 
are localized simultaneously on second- and first-layer 
sites, for example, on Ti(32) and Ti(6) or on Ti(21) and 
Ti(15). Furthermore, a configuration has been sampled 
where the excess charge is shared between the second 
and third layer involving Ti(30) and Ti(40). The av- 
erage lifetime for a specific charge-localized topology is 
roughly 0.3 ps at 1000 K, which underscores the pro- 
nounced dynamical nature of the phenomenon. At 700 K, 
this timescale increases to about 0.4 ps, from which naive 
Arrhenius extrapolation yields a formal activation energy 
of approximately 0.07 eV and lifetimes of the order of ps 
and ns to /is at ambient and liquid nitrogen temperatures, 



respectively. This implies that dynamical averaging up 
to quenched disorder at liquid helium conditions might 
be operational in surface science experiments. 

Next we consider the complementary scenario where 
the defect charge is localized initially in both the first 
and second layers at Ti(ll) and Ti(29) (see Fig. GJb)). 
After only 0.1-0.2 ps, the charge at Ti(ll) transfers from 
the Ob surface row to the adjacent surface row of Tis c 
atoms where it is delocalized over two nearest-neighbor 
sites Ti(14) and Ti(15). It remains shared between these 
two sites for about 0.5 ps before it jumps to the nearest- 
neighbor second-layer atom, Ti(31). Here, charge trans- 
fer between the first and second layers is found to be me- 
diated by a charge derealization process involving first- 
layer Ti sites in the Ti5 C row. Once the charge reaches 
the second layer, during the next ~ 4 ps it visits the same 
Ti sites as it had earlier explored. In one of these events 
a flipping of one of the two unpaired electrons between 
an a- and (3-d orbital occurs along the trajectory. 

The dynamical behavior can be cast into distribution 
functions using localization histograms, which is illus- 
trated in Fig. [21(c) for the dynamics in (a). In both set- 
ups the excess charge is found to be extremly mobile and 
to eventually visit all available Ti sites, although it dis- 
plays a strong preference for populating sites in the sec- 
ond subsurface layer under the Tis c rows (roughly 70 %). 
This is followed by population of surface states 20 %), 
and localization in the third subsurface layer (« 10 %). 
We note that in all observed charge transfer events there 
was always a mediating nearest-neighbor Ti site tran- 
siently involved in the process. This dynamical scenario 
is in substantial agreement with the viewpoint suggested 



based on recent experiments 

In order to assess the energetics of the various local- 
ization topologies, we have characterized representative 
structures in detail. Six additional excess charge topolo- 
gies were obtained by quenching a large set of sampled 
configurations to T = K using Car-Parrinello anneal- 
ing [24[ followed by standard optimization (see Table H]). 
Localization of one excess electron at Ti(ll) induces a 
significant elongation of the Ti 3+ fl l)-0 distance (see Ta- 
ble Ij), which is as expected 15, 0, [l9[. The relative 
energies confirm that the most stable sites for charge 
localization are those where both electrons are in the 
second subsurface layer under Tis c rows - in line with 
earlier findings [2l[ and consistent with our localization 
histograms. About 0.2-0.3 eV higher in energy are the 
topologies where the excess charge is shared between sur- 
face Ti 5c atoms and second subsurface layer sites below 
Ti5 C rows. Topologies resulting in excess charge local- 
ization on both second and third subsurface layer sites 
under Ti5 C are less stable by about 0.3-0.4 eV. In the 
least stable configuration, (by 0.7-0.8 eV) one electron is 
at an under-coordinated surface site, Ti(ll) at the Ob 
vacancy, whereas the other electron has transferred to 
the second subsurface layer under a Tis c row at Ti(29). 



In one case, starting from a configuration in which the 
charge is delocalized over two surface atoms, Ti(14) and 
Ti(15), and one second layer site, Ti(29), we end up with 
a solution in which one electron is trapped at Ti(29) un- 
der the Ti5 C row and the other electron is localized on the 
Ti site directly below the Ob vacancy. This configuration, 
which has also been recently reported [l6| , is found to be 
0.22 eV higher in energy than our lowest energy configu- 
ration. The concomitant static distortions obtained from 
the nearest-neighbor Ti 3+ -0 distances of the optimized 
structures (b)-(d) in Fig. [T] relative to the non-relaxed 
reference (a) are typically 0.05-0.1 A. However, already 
at 300 K the thermal fluctuations are found to override 
these static distortions by inducing larger root-mean- 
square-deviations of the same relative distances. 

Combining our static and dynamic calculations demon- 
strates that the excess charge dynamics is intimately 
tied to the existence of many close-lying minima on the 
PES which can be populated at finite temperatures [12 1. 
Charge localization pattern from static optimizations 
therefore depend heavily on the initial configuration due 
to trapping in local energy minima. This might con- 
tribute to the variety of differing results reported in the 
literature. In this situation, AIMD helps greatly by not 
only providing a dynamic perspective, but also by being 
able to explore vastly different localization topologies. 

Finally, the total and projected electronic density of 
states (DOS) of the distinctly different charge localiza- 
tion scenarios are analyzed in Fig. Efd). In all cases, 
the O vacancy leads to the appearance of distinct, filled 
states in the TiO2(110) band gap, whose projection onto 
two specific Ti sites - although they are different for each 
of the three scenarios - fully accounts for the entire peak 
integral area. Together with the almost instantaneous 
changes in localization patterns embodied in Fig. [2ja) 
and (b), this indicates that the two excess electrons are 
essentially always trapped at two well-identified, specific 
Ti sites. These sites, however, are not unique for a given 
defect, but rather they interchange dynamically. 

In conclusion, finite temperature PBE + U simulations 
paint a complex picture of the dynamics of electrons in 
defects at the reduced titania TiO2(110) surface. The ex- 
cess charge, being trapped at specific Ti sites, migrates 
easily by phonon-assisted (thermally activated) hopping 
to other Ti sites, thus exploring significantly different 
electronic structure topologies. Due to averaging effects, 
this leads to an "effective derealization" , which is pref- 
erentially found on second-layer subsurface Ti atoms un- 
derneath Ti5 C rows. Although much less frequently, ex- 
cess charge also "visits" first layer and third layer subsur- 
face sites. We expect that this defect-induced complex 
charge (de-) localization dynamics scenario is not only 
fundamental to titania, but of broad significance to ox- 
ide materials in general. 
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